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Abstract 

We analyze the dynamics of a deterministic model of inhibitory neuronal networks 
proving that the discontinuities of the Poincare map produce a never empty chaotic 
jy^ • set, while its continuity pieces produce stable orbits. We classify the systems in three 

^"*^ | types: the almost everywhere (a. e.) chaotic, the a. e. stable, and the combined systems. 
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The a.e. stable are periodic and chaos appears as bifurcations. We prove that a.e. 
stable systems exhibit limit cycles, attracting a.e. the orbits. 

Keywords: Chaotic sets, limit cycles, piecewise continuous dynamics, neuronal net- 
works. 

INTRODUCTION 

We obtain rigorous mathematical results on the dynamics of a deterministic abstract 
discontinuous dynamical system, in a finite but large dimensional phase space. It comes 
from a non linear model of inhibitory neuronal networks, without delays, composed by 
equally or different 2 < n < +oo pacemaker neurons, evolving according to an autonomous 
■ differential equation in the inter-spike interval times, and interacting among them by 

synaptical instantaneous currents in the spiking instants. 

A vectorial autonomous differential equation governs the increasing potentials of the n 
neurons as a function on time t, only during the inter-spike regime. On the other hand, 
^ . the spiking regime holds when at least one neuron, say i, reaches a threshold level and 

gives a spike. Due to the synaptic connections among the neurons, this spike produces 
sudden changes in the potentials of the other neurons j / i and resets the potential of the 
neuron i. 

The synapsis is assumed to be inhibitory, i.e. phase redeeming, meaning that the potentials 
of the receiving neurons suffer negative changes in the spiking times. That is why the 
inhibitory synapsis is modeled by a matrix {Hij}^, 1 < i,j < n of negative numbers 
Hij < that represent the instantaneous discontinuity jumps in the potentials of the 
neurons j ^ i, produced by a spiking of the neuron i. 

The autonomous differential equation, verifying some very wide assumptions, governs the 
system during the inter-spike intervals of time. It leads to a Poincare map (Sotomayor 
[1979]), which is contractive, as we prove in Theorem 0] (see also Budelli et aL[1996]). 
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Finally, when synaptic inhibitory coupling is added, it produces asymptotic cyclic at- 
tractors in the phase space, whose existence we prove in Theorem [U]) . Nevertheless, the 
discontinuities of the synapsis also generate chaotic orbits, whose existence we prove in 
Theorem [71 While the limit cycles exhibit mostly asynchronous spikes, but in periodic 
patterns, they could seem irregular to an experimental observer, if their periods are very 
large, and be virtually mixed with the real chaotic orbits. 

The model we study includes the networks of coupled leaky integrators and the relaxation 
oscillators. The mathematical approach, under particular hypothesis, on the integrate and 
fire and the relaxation oscillator models of many neurons was early analyzed by Mirollo 
& Strogatz [1990] for homogeneous networks of excitatory neurons. Later Abbott & van 
Vreeswijk [1993] generalized the results for non homogeneous systems, also with delay 
times and excitatory interactions. For strongly coupled networks of inhibitory or excitatory 
neurons Bessloff & Coombes [2000] find synchronous and asynchronous behavior. 
In the abstract dynamical system we define the complementary sets C and S, of chaotic 
and stable future orbits respectively (Section 3). 

The qualitative rigorous analysis, discrete-sizing the dynamics by means of a Poincare 
map (see Sotomayor [1979]), was applied to homogeneous networks by Mirollo & Strogatz 
[1990]. Also Budelli et al. [1991] and Catsigeras & Budelli [1992] studied the asymp- 
totic future behavior of inhomogeneous two neurons networks, with a Poincare map in a 
co dimension one section. 

On the other hand discrete-sizings on time, are different mathematical approaches and 
can fit better with computer experiments, in which both time and space are discrete- 
sized. Cessac [2008] and Cessac & Vieville [2008] obtain similar results to those in this 
paper, using different mathematical models. They also study discrete discontinuous neural 
networks. Their model, in spite of being similar in the evolution given by equation ([I]), is 
intrinsically different to consider the spike instants as predetermined by the observer. 
In this paper we prove that the Poincare map is discontinuous due to the inhibitory 
synapsis. Although it is piecewise contractive, we prove that its discontinuities play the 
role of chaos generators: due to them the chaotic set C is never empty. 
The abstract dynamical systems of this model are classified in three types according to their 
attractors a.e.. Calling m to the Lebesgue probability measure, we define the a.e. chaotic 
systems (for which C has full measure: m{C) = 1, m{S) = 0); the a.e. stable systems (for 
which m(C) = 0, m(S) = 1); and the combined systems (for which m(C) > 0, m(S) > 0). 
Generically, under some additional hypothesis, the systems are a.e. stable, and the chaotic 
set appears as a bifurcation among different stable systems. This kind of bifurcations were 
studied by Catsigeras, Rovella & Budelli [2008], proving that they exhibit a chaotic Cantor 
set attractor. 

In Theorem [9] of this paper we prove that the a.e. stable systems exhibit limit cycles 
attracting Lebesgue almost all the orbits, and due to them they are a.e. periodic. Nev- 
ertheless the period of the attracting limit cycles, and the number of them, are mostly 
determined by the relation among the inhibitory synaptic interactions. They can a priori 
be very large, and have few relation with the intrinsic individual periods of the integrate 
and fire neurons of the network. This last result fits with the numerical experiments in 
networks of n different integrate and fire neurons. Postnova et al. [2007] analyzed the 
dynamics of computer simulated large networks, and obtained that the system can be 
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driven through different synchronization states, but they are significantly different from 
the original periodic behavior of the individual cells. 

1. DEFINITIONS AND STATEMENTS. 

The phase space is a bounded cube Q = [— 6,6] n of the vectorial space M n , whose points 
are the n-uples V = (V%, . . . , V n ) with — 9 < Vi < 9, V 1 < i < n. We denote 

I = {l,2,...,n} 

to the set of neurons of the network. The variable Vi is the potential states of the i— th. 
neuron of the network, and it evolves on time t. 

In this model the synapses and the spikes of the neurons, can be produced at any instant 
tf. when the potential Vi of at least one of the neurons (say the i-th neuron), reaches 9 
(the fixed threshold level), ij is not previously specified: it is the solution of the implicit 
function = 9. We denote as to the i-th. component of the vectorial flux solution $ 
of an autonomous differential equation 

dV/dt = F(V), V7eQcK n , F : Q ^ R n , F G C 1 (1) 

This vectorial differential equation determines the potential Vi of each neuron i as a func- 
tion on time t, only during the inter-spike regime. The solution <&*(V°) depends on the 
initial state V° = (Vf, . . . , V°). We restrict to the case in which the system is composed by 
n independent differential equations (i.e. Fi depends only of Vi). Nevertheless the open 
conclusions we obtain (for instance the result in Theorem U|) holds even if the neurons 
slightly interact during the inter-spike intervals (i.e. dFi/dVj ~ if i ^ j). We assume 
that 

Fi > 0, dFi/dV < Vie/, dFi/dVj = V;?' ^ i (2) 

and Fi and dFi/dVi are bounded away from zero. The dynamics given by ([I]) and ([2]) 
is that of general pacemaker neurons, with strictly increasing potentials in dissipative 
regime, providing negative concavity of the free evolution on time, of the potential of each 
neuron in the inter-spike time intervals. This property leads to a contractive Poincare 
map (Theorem [4} and produce, when synaptic inhibitory coupling is added, asymptotic 
cyclic attractors in the phase space (see Theorem [9|), but also chaotic orbits (see Theorem 
[7]). The limit cycles exhibit mostly asynchronous spikes, but in periodic patterns. 
The model verifying ([1]) and ([2]) includes the networks of coupled leaky integrators and 
the relaxation oscillators, for instance if Fi = — ji(Vi — /%), with constant < 7i, ft. 
From an initial state V°, the next spiking instant is 

t(V°) = min U(V°) (3) 

l<j<n 

where U(V°) is the solution of the implicit function 

$*(y°) = 9 ^ t = ti{V°) (4) 

At instant t maybe more than one neuron reaches the threshold level 9 simultaneously 
(although this occurs with zero Lebesgue probability in the initial state, if the system is 
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inhibitory). We call J(V ) to that set of neurons that spike first, from the initial state 
V°, i.e. 

j(y°) = {ieI:t i (V°) = t(V )} (5) 

Almost all initial states are such that J(V°) = {%} for a single neuron i G /. The set J of 
neurons gives spikes at time t, i.e. their potentials reset to zero, and from the new initial 
state W° it restart the evolution according to the differential equation ([I]). Precisely 

Vi{T) = 9 => Vi(t + ) = Wi(0) = Vie J (6) 

When spiking, each of the neurons i £ J, which are supposed to be inhibitory, produces a 
sudden synaptic current through its connections to the other neurons j €" J. This synapsis 
produces a sudden change, of amplitude —Hij < 0, in the potential Vj. Precisely: 

Vj(F) < o, Vj(t + ) = wm = VjCn - H u V ? J (7) 

2. THE POINCARE SECTION AND THE FIRST RETURN MAP. 

In this section the dynamics of the first return map p : B i-> B is adequately defined in a 
Poincare section B (see the definition of Poincare section and map in Sotomayor [1979]). 
The section B will be diffeomorphic to the union of a finite number of n — 1 dimensional 
balls transverse to the flux. 

In the compact phase space Q = {V G M" : —9 < Vj < 9 V i € J} we take the following 
n — 1-dimensional section: 

n 

B=[jB k where B k = {V € Q : V k = 0} (8) 

k=l 

Its topology is that induced by B k C M n_1 . 
We assert that B is transversal to the flux $ . 

Proof: It is enough to prove that is transversal to the flux for all k = 1, . . . , n. This 
last assertion is deduced from the property d<&\/dt = > 0, due to ©, as follows: the 

n components of the vector d&/dt, which is tangent to the flux, are positive, in particular 
its k— th. component. On the other hand, the manifold B k is n — 1 dimensional, defined 
by the equation V k = 0. Therefore its tangent subspace S k ~ M n_1 is formed by all the 
vectors in R n such that have null their k-th. component. Then we deduce: 
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denotes the subspace generated by the vector U G R n . □ 

Due to the definition of the mathematical spike in equation ([6|), from any initial state 
V° £ Q the system arrives to the Poincare section B in a finite time t(V°) given by 
equations ([3]) and dH). So in particular for V° G B, we have defined the first return 
Poincare map p: 

p:B^B, p(V°) = W° = V I+ ^ (9) 
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The sequence of inter-spike intervals (ISI) is given by the evolution of the system through 
the iterates pP,p = 1,2,... of the return map p. The ISI sequence is 

t(V°),t(p(V )),t(p 2 (V°)),...,t(ff(V )),.... 

If the orbits of the discrete dynamical system by iteration of p are attracted to limit cycles 
or not, the ISI sequence will be finally periodic or not. If the dynamics of p exhibits chaotic 
attractors, then the ISI sequence also. Even in the periodic case, its period depends of 
the map p, (and also of the initial state if there were many limit cycles). But p depends 
strongly of the matrix of synaptic interactions Hij. Therefore, the network composed by 
the inhibitory coupling of n oscillators, may have a dynamics which widely differs from 
the behavior of each isolated neuron. 

We iterate the Poincare map p(V ) = W°, and after each iterate we reset the time to 
consider the solution $*(W°) of the equation (pQ), from the new initial state W° G B. To 
simplify the notation we will omit the supraindex °: from now U,V,W denote points in 
the Poincare section. 

From equalities ©, ([7]) and ([9]) we have the following formula for the Poincare map: 

P 0)= *5 (y) a/)-£ ie j ( v)i^ vj^j(y) (io) 

Pi {V) =0 Vi G J(V) 

Denote p{I) to the family of all non empty parts of the set / of neurons. Roughly speaking 
p(I) the collection of all possible "words" of different neurons, of any length > 1, < n, 
without considering the order. For each J G p(I) (say for instance, J = {1,2}), define: 

Bj = {V G B : J(V) = J} 

In the example J = {1,2}, the set -B{i,2} C B is composed by all the initial states in 
the Poincare section B for which neurons 1 and 2, and only them, will arrive first and 
simultaneously to the threshold level. 

Remark 1 Consider the difference, not only in notation but in significance, among the 
sets B^j and Bi. 

On one hand, B^y, for a fixed index i = 1,2, ... ,n, is the set of all the initial states in 
the Poincare section B, from which neuron i, and only neuron i, will arrive first to the 
threshold level, after the system's evolution in the inter-spike regime. Then BuyHBu-i = 
if i + j. Also Uf =l B {l} C^B = Uj €nP{I) Bj. 

On the other hand, according to its definition in ([8]), B^, for a fixed index k = 1, 2, . . . , n, 
is the set of all the states in the Poincare section B, in which the potential level V*. of the 
neuron k is the reset level zero. We have B^V\Bh = {V G B : Vh = Vk = 0} ^ V /i / k, 
and B = U^^B^. Also, due to the definition of the Poincare map p in the equality after 
CED, we deduce that p{B^) C Bi, p(B {i j y) cB ; n Bj. 

Remark 2 Due to formulae (I10|) the Poincare map p is continuous in each piece Bj, in 
particular in B^y for all i G I. The set Bj, for each J G p(I), is called a continuity piece 
of p. 

From (|10p . due to the transversality of the flux <I>* to Bi, p transforms homeomorphically 
B^y n Bk onto its image in Bi. So it is an open map. Nevertheless it is not globally 
injective because each point may have different pre-images in Bk and Bh for h ^ k. 
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The partition V and the set dV of discontinuities. 

By construction, for different subsets J 1 ^ J 2 G p(I)- Bj 1 f\Bj 2 = 0. Any initial state 
V in B belongs to one and only one set Bj (exactly to that Bj such that J(V) = J). 
Therefore: 

The collection {Bj} jg^m form a partition V of the Poincare section B into the continuity 
pieces of the return map p. 

© and @ implies that the set of V G B such that t(V) = U(V) < tj(V) Vj ^ i is open. 
In other words it is open the set B^ (i.e. the set of initial states from where the neuron 
i, and only the neuron i, arrives first to the threshold level). On the other hand Bj has 
empty interior if # J > 2. 

We consider only those not dead neurons (i.e. they do not remain under the threshold 
level forever) and rename as / to that set of neurons. If there is only one, then there is not 
interactive network to study. So let us suppose that > 2 and then {Bi ^ : i G 1} 
is a collection of at least two open and non empty sets. As the space B is connected, we 
conclude that the topological frontiers of the pieces of the partition V is not empty. Define: 

n n 

dV = B\[jB {i} = |J Bj = {JdB {l} . 

*=1 #J>2 *=1 

Topological properties of the Poincare map. 

We use the following norm to compute the distances in each Bk C B, for k € /: 

HV-Wl^maxl^-Wil VV,WeB k (11) 

Assume the following generic hypothesis in the parameters of the system: 

Hij^e Vi^jel 

Define the expansivity constant: 

min,-/i \8 — Hi A 
a = 11 > (12) 



Lemma 3IfV£ dV then the Poincare map p is discontinuous in V and the discontinuity 
jumps in V are larger than 3a. Precisely, there exists a sequence of points U m — > V such 
that lim m ^ +00 \\p(U m ) - p(V)\\ > 3a 

Proof: If V G dV then V G Bj for some J G p(I), #«/ > 2. There exist i ^ j in J and 
therefore Bj C dB^. By the definition of topological frontier, there exists a sequence of 
points U m -> V such that U m G B {i} . 

Pj (U m ) = <5>} Um) (U m ) - -^ m ^ + oc ®} V) {V) - Hij 

The last limit is computed recalling that the solution of the differential equation ([T|) 
depends continuously on the initial state V, and the ISI time t(V), defined as the minimum 
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in equality ([3]) of the solutions of the implicit equations (jl]) , also depends continuously on 
V, due to the Implicit Function Theorem. (Rey Pastor & col. [1968]). 
But j G J, V G -Bj ^ j £ ^(^)) or m other terms, V is an initial state for which 
the neuron j arrives to the threshold level. Then $*• (F) = 6, and by equalities (fTOj) . 
j e J(V) ^(V) = 0. Then: 

limpet/™) = 0-%, p j (V) = 
| lim (C/ m ) - Pj(V)\ = \9- Hij | > 4a > 3a 

By the definition of limit, for all m large enough we obtain \\pj(U m ) — Pj(V)\\ > 3a. □ 
The Lemma [3] explains why the discontinuities of p act as saddle type orbits in continuous 
systems, to produce expansion and chaos. The points of dV act on their neighborhoods 
with an infinite rate of expansion in some directions, given by instantaneous expansivity 
larger than an uniform number 3a > 0. But they also have directions of contraction, as 
they are in the frontier of the open pieces of continuity of p, where the following 
result holds: 

Theorem 4 The Poincare map p : B i— > B is uniformly contractive in each of its conti- 
nuity pieces B^ n p(B) . Precisely, there exists < A < 1 and a distance dist c such that, 
ifV,U G -B{j} n p(B) for some i £ I, then 

dist c (p(V),p{U)) < Xdist c (V,U) 

Besides, there exists K > such that for all V,U G B: 

— \\V - U\\ < dist c (V, U) < K\\V — U\\ 

We note that the distance dist c is not necessarily the norm || - 1| defined in (jlip . Neverthe- 
less, they are equivalent in M ra_1 . For simplicity we will omit the not zero constant factors 
1/K and K in the forward computations along this work, and use |[ • |[ as if it were dist c . 
This is not a restriction because, as a consequence of Theorem 01 for any norm there ex- 
ists an iterate po such that (f° is contractive in its continuities pieces: (^(V) — fP(U)\\ < 
Kdist c {pP{V),pP(U)) < K\Pdist c (V,U) < K 2 \p\\V - U\\ < (1/2)||V - U\\ for all p large 
enough. 

Proof: Define < eo = min^j H^j- Apply equality (fTUj) and use, at the end of the 
following formulae, the Lagrange Mean Value Theorem (Rey Pastor & col. [1968]), to 
compute the difference of $*(V), for different values of the time t, as the increment At 
multiplied by the derivative in an intermediate time T m : 

Ve P (B {k} ) V k = 0, Vj<e-e VjeI 

V G B k = {V G B : V k = 0} 

VGB {1} $ V \V) = 6, $(y) = Vi 

e <e-v i = t{v)d®\(y)/dt\ t=Tm 

eo<t{V)F^ Tm {V)) 
=> t(V) > e /(maxi<i< n maxy e q Fi(V)) = t > (13) 
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By the Tubular Flux Theorem (Sotomayor [1979]), there exists a C 1 bounded diffeomor- 
phism which is a spatial change of variables £ : V i— y V from Q C W 1 onto Q C M n , such 
that £\b = id and the solutions of the differential equation (P) in Q verify dV/dt = a in 
Q, where a £ W 1 is a constant vector with positive components. It verifies: (V)) = 
COO + a-i, cZ£F(V) = a VKgQ. 

Define in R n the orthogonal projection ir onto the (n — l)-dimensional subspace a\V\ + 
«2^2 + • • • + a n V n = 0. The flux of the differential equation ([1]), after the change £ of 
variables in the space, is ortogonal to that subspace, and is transversal to B^. Consider 
any real function g: 

VV,V + dV£R n : 
Tr(dV) = Tt(dV + ag(V, dV)). 
VV, V + dV, U G B fcl define 
dist c (T/,y + dV) = \\w(d£vdV)\\ 

dist c (V, CO = / Ik ■ dev +t( u-v) • (t/ - V) || tft (14) 
■/ o 

As £ is a C 1 diffeomorphism, its derivative and the derivatives of its inverse, are bounded 
in the compact set B, and so the distance dist c defined above verifies the last thesis of this 
Theorem. It is left to prove that p is contractive with this distance in B^ D B^. 
Let us apply p to V, V + dV £ B^y n B k . We use the equality with J(F) = {i}, in 
which for convenience we agree to define Ha = 6. We shall use the derivation formula of 
the flux of the differential equation respect to its initial state: 

d^ j {V j )/dV j = exp (ffaFj/dVj) (^(Vj)) dsj) 
Define: — 7 = maxi<j maxy £ Q dFj(Vj)/dVj < 

In what follows i is fixed. It is the value of the index of the continuity piece B^y given 
in the hypothesis of this lemma. On the other hand, j = 1, 2, . . . , i, . . . ,n is the index of 
the general component dQydt of the tangent vector of the flux $*, and of the general 
component Pj(V) of the Poincare map whose derivative we are computing. We must 
compute all their components, so we must include the case in which j = i. 
The formula (fit)]) and the inequality (fT3|) lead to: 

p(V)-p(V + dV) = dp(V)-dV = [(dpj(V)/dVj)dVj + {dp j {V)/dV{)dV i \ l ^ n dp j {V)/dV j 

= (^(^)/^)| t= j (y) =exp n* iV \dF j /dV j )($ s j (V j ))ds)j <e-*>d Pj (V)/&Vi = 

= (dmVj)/dt)\ t=m ■ (dt(Vi)/<W) = g(V) ■ F^XV)), 

where g(V) = dt(Vi) / dVi is the real function obtained deriving respect to Vi the implicit 

equation 6 = \Vi). Call ej to the j— th. vector of the canonic base in W 1 and join all 
the results above: 

7T • di p{y) ■ ( P (v) - P {v + dv)) = vr ■ di piy) ■ d P (v) ■ dv = 
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* ■ d Zp(V) ■ YtVPityyaVi) ■ dVjgj )+ n ■ di p{v) • (g(V) ■ F(^ V \V)) 



J =1 



J =1 



* ■ d Uv) ■ [ Y^idPjOO/dVj) ■ dV 3 e} ) + n(g(V) ■ 

= vr • d6p(V) • (j2(d Pj (V)/dV 3 ) ■ dV 3 ejj (15) 

Now we define the number < A = e~ 7 * < 1 and observe from the computations above 
that: 

< dp{V)/dVj < e~ 7 * = A < 1. 



Applying the definition of the differential distance dist c in (|14p , and the equality (|15p , we 
obtain: 

dist c (p(V),p(V + dV)) = \Wdt p(v) dpv -dV)\\ < 
A \\7rdC ■ dV\\ = A dist c (y, V + dV) = A ||7r(d£ydF) || 

Integrating by formula (|14|) we conclude: 

dist c (p(V),p{U)) < Adist c (V,*7) □ 

Measure properties of the Poincare map. 

Let m be the (n — l)-dimensional Lebesgue probability measure in the Poincare section 
B. Let us prove that 

ro(^)=m(U #J > 2 Bj)=0. 
In fact, each Bj, if # J > 2, is the finite union of C 2 n. — J)— dimensional manifolds in 
B, obtained when, for two or more different values of i £ /, the respective solutions ti(V) 
of the implicit equations (j4|) , are equal. 

Our aim is to study the attractors. The weakest condition required to a set A C B to 
be an attractor is that its basin of attraction has positive Lebesgue measure. So, we may 
take out the points of the measure zero set dV. 

The set B' of the points with infinite itinerary. 

Define 

B' = {V £ B : pi{V) £dV Vj > 0} C B \ dV 

B {i} B {J} = if M 3 and UieJ B {i) = B \dV => 

5' = n p + = °°o ^ P (^ \ flP), B\B' = U+ = °° 

The following assertion characterizes the set B': 

V G B' 4$ 3 a unique sequence {i p } p >o such that pPiy) € B^ ip y \/p > 0. 
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This sequence is called the itinerary of V , and it is the infinite sequence of neurons that 
will spike (in different times), in the order they reach the threshold potential, from the 
initial state V. 

p is continuous in each Buy Therefore, for all p > and for all V £ B', the iterate map 
pP is continuous in V. As B' is the numerable intersection of open and dense sets, it is 
dense. 

Theorem 5 The set B' with infinite itinerary has full Lebesgue measure in B, i.e: m(B\ 
B') = 0. 

Proof: To prove that m(p^ 1 (dV)) = 0, once it is known that m{dV) = 0, use the same 
argument as to prove by induction that for all p > 1 : m(p~ p (dV)) = 0, once it is known 
for p — 1. 

In fact, we apply Liouville Formula (Sotomayor [1979]) to compute the following Lebesgue 
integral: 

= m(dV) > m{p{p~ 1 (dV))) = 

|det(-Dp)(F)| dm(V) 

Vep-^i&P) 

and conclude that m(( / 9~ 1 {dV))C\{V £ B; [det (Dp) | > 0}) = 0. To prove that m(p- 1 (dT>)) ■ 
it is enough to prove that the Jacobian det(Dp) ^ a.e. in B. 

To apply the Liouville Formula we shall first prove that p, given by formulas (|10p . is 
differentiable a.e. in the integration set. The first technical problem arises because p is 
not differentiable in all the points: it is neither in the points of discontinuity dV nor in 
the points of \Jh^kel(^h H B^). In fact, in these last set the Poincare section itself fails to 
be a local differentiable manifold. The set of those exceptional points has zero Lebesgue 
measure because they are contained in (n — 2)-dimensional subspaces. 
The second technical difficulty is to check that det Dp ^ a.e. The derivative of the 
Poincare map can be computed directly from formulas (|10p . in each continuity piece B^ 
intersected with each B]~ \ [J^^B^. In those sets B is a local differentiable manifold. 
Computations in (??) lead to 

det Dp(V) = Yij^i j^k (9pj(V)/dVj) ■ (dpk(V)/dVi) All these factors are computed in 
equalities (??) and are not zero due to hypothesis ([2]), so det Dp ^ a.e. □ 

3. CHAOTIC AND STABLE SETS. 

We will divide the Poincare section B in two complementary sets S, formed by stable 
orbits, and C = B\S, the chaotic ones, according to Definition [U The set S is formed by 
stable future orbits under any uniform sufficiently small perturbation that can be added at 
any step of the iteration of p. In the set C = B\S there are arbitrarily small perturbations 
that, if added in some step of the iteration of p, drastically change the future orbits and 
their asymptotic behavior. To have a criteria of chaos, what we call drastic changes in the 
phase state, we consider the expansivity constant a > defined in (I12|) . 

Definition 6 Stable and chaotic sets. V € B is stable if there exists 5 > such that 
Vp > 0, \/W G B, if \\pP(V) - W\\ < S then \\p k (p p {V)) - p k {W)\\ < aVfc > 1. 
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V G B is chaotic if it is not stable. The opposite of the definition of stable point holds as 
follows: V G B is chaotic if and only if for all 5 > there exists p > 0, and there exists 
W eB such that \\pP(V) - W\\ < 5 and \\p k (p p (V)) - p k (W)\\ > a for some k > 1. 

5 is i/ie set o/ t/ie stable orbits and C = B\S is the set of all the chaotic orbits. 

It is immediate from the definitions above that S is forward invariant: p(S) C S, and 

thus, its complement C is backward invariant: p _1 (C) C C . 

Given 5 > fixed, we define the uniform stable set S$ C S as the set of V G B such that 
Vp > 0, W G 5, if ||pP(V) - W|| < <5 (with 5 previously fixed), then \\p k (p p (V)) - 
p k (W)\\ < a\/k > 1. From the definitions above observe that p(Ss) C 5<s and £/zai any 
point in S is in Sg for some 5 > 0. Therefore, to study the dynamics of the stable points, 
it is enough to study the dynamics in the sets S$. 

Theorem 7 The set C of chaotic points is never empty. Precisely: C D B\B' D dV ^ 0, 
where B' is the set of points in the Poincare section B with infinite itinerary, and dV is 
the set of discontinuities of p. 

Even more, if V G dV then for all 5 > there exists W G B such that \\V — W\\ < 
5, \\p(V) — p(W)\\ > 3a, where a is the expansivity constant of the system. 

Proof: From the note at the end of Remark [2] and from the definition of the set B' of the 
points with infinite itinerary, we get ^ dV C B \ B'. 

First, let V G dV. We shall prove that V G C. From Lemma [3] there exists sequence of 
points U m -> V such that lim \\p(U m ) - p(V)\\ > 3a. Given 5 > for all m > 1 large 
enough the points U m — > V verify \\U m — V\\ < S. Therefore, taking W = U m the point 

V verifies the definition of chaotic point with p = and k = 1. 

Second, take V G B \ B' = \J p>0 p- p (dV). We shall prove that V G C. We know that 
for some p > p^iV) G dV. In the first step we proved ^(V) G C, so V G p~ p (C) C C, 
because C is backward invariant. □ 

Classification of systems. 

► a. e. stable systems: m(S) = 1 and m(C) = 0. In Theorem [9] we prove that the a.e 
the limit set is composed only by periodic sinks (limit cycles). Nevertheless the set C of 
chaotic points in not empty as proved in Theorem [3 

► Chaotic systems: m{C) = 1 and m(S) = 0. By definition the limit set A of C is a 
chaotic attractor. In Catsigeras & al. [2008] it is proved that A is a Cantor set attractor. 

► Partially chaotic and partially periodic systems: m(C) > and m(S) > 0. We will 
show that the points in S are attracted to limit cycles. 

Remark: We did not find C 1 examples in our model of systems for which m(C) > 
although due to Theorem [7] the set C of chaotic points is never empty. Nevertheless it is 
possible (but not immediate) to construct piecewise continuous C° maps pinan — 1>2 
compact ball B, that are uniformly contractive in each of their continuity pieces and such 
that C = B, 5 = 0. 
Some known results. 

Suppose allowed any small perturbation of the system structure, such that, instead of 
having the differential vectorial equation ([T]), we have dV/dt = G(V) with ||G — < e 
for some e > sufficiently small. By continuity in the C 1 topology of the functional 
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space, the first two assumptions ([2]) are also verified by G and the third one transforms in 
dGi/dVj « Vj^i. 

Besides, instead of restricting to a constant matrix Hq = [Hjj]jj e j to describe the synaptic 
interactions, allow a matrix —H(V) of negative numbers —H(-) < to be a C 1 function 
such that — -frolic 1 < e - 

The new Poincare map p will move respect to the old one. The new set of discontinuities 
dV of the Poincare map will be (with the Hausdorff distance) near the old one, because 
its points are those defined by implicit functions, whose equations have C 1 dependence on 
G,H £ C 1 . So they are near the old ones if \\G — F\\ c i and \\H — -ffolle 1 are small enough. 
We conclude that, if e > is sufficiently small, the new Poincare map will be C 1 near the 
old one, in each of its continuity pieces. 

In this scenario, the following results are known, leading to the generality of the property 
m(S) = l 5 m(C) = 0. 

► If provided the additional open hypothesis of separation (pairwise disjointness of the 
closure of the images of the different pieces of continuity for some iterate of p), then 
Catsigeras & al.[2008] proved that generically (open and dense in the C° topology), there 
exist a finite number of limit cycles that attracts all the orbits of B' . We conclude that 
S = B', and thus m(S) = 1. 

► If the contractive map p is affine in each of its continuity pieces, then Cessac [2008] 
proved that for generic values of the real parameters, not only in the topological sense 
(open and dense family of systems), but also in the Lebesgue measure sense (in the space 
of finite number of the real affinity parameters), there exists at least one and up to a finite 
number of limit cycles, which attract all the orbits of B'. We conclude that S = B' and 
m(S) = 1. 

► If the contractive map p is not affine, but the matrix of interactions H is constant, 
Cessac & Vieville [2008] proved that either there exists limit cycles attracting the points 
of B' or the dynamics has positive entropy. 

Definition 8 Omega limit and limit cycles. 

For any V G B, its omega-limit co(V) G B is the set of limit points of the future orbit of 
V. Precisely: 

u(V) = {W eB : 3 Pj -> +oo such that pP* (V) -> W} 

As the phase space B is compact, the omega limit set of any point is not empty. The 
Poincare map p is not continuous, so oj(V) is not necessarily forward invariant. Neverthe- 
less the omega limit set is the same for all the points in the same orbit, as it is easy to 
check. 

A set L G B is a limit cycle (also called a periodic sink) if it is a single periodic orbit and 
its basin attraction B(L) contains an open neighborhood of L, being: 

L = {V,p(V),...,p r - 1 (V)}, P r (V) = V r>\ 
B{L) = {W G B : u(W) = L} 

Theorem 9 If for some 5 > the set Ss of uniform stable points is not empty then, there 
exists a finite number N of limit cycles L\, . . . such that the union UfcLi B(Lk) °f 
their basins of attractions includes Ss- 
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Remark: Let us show the main consequence of this theorem: if the set S of stable points 
has positive measure, it is not empty and, because all its points are included in uniformly 
stable sets Ss, they all are attracted to limit cycles. By definition of limit cycles, the 
basins of attraction of different limit cycles are pairwise disjoint containing open sets. The 
topology in B has numerable basis, so we conclude that there are at most a numerable 
quantity of limit cycles attracting all the orbits of S. If besides m(S) = 1, then a.e. point 
in B is attracted to a limit cycle. 

We will prove Theorem [9] at the end of the next Section. 
4. THE ATOMIZATION OF THE SPACE. 

Let us consider a system for which the stable set S is not empty. All the results in this 
Section will hold under this hypothesis. 

As observed at the end of Definition [6) S ^ => Sg / for some 5 > 0, and thus also 
for all sufficiently small 5 > 0. Fix < 5 < a such that Ss ^ 0, where a is the expansivity 
constant. We recall that Ss C B' from Theorem [71 and p(Ss) C 5,5. Therefore the image 
^(Ss) is disjoint with dV for all p > 0. So it is partitioned in n disjoint pieces (may be 
some of them are empty) when intersecting Ss with {B^y}i e i. In each of this pieces p is 
continuous, and p(Ss) = (Jig/ p(Ss) H ^{i}- 

Definition 10 Atoms of generation p. Given i\ E / we call atom of generation 1 to 
the following set: 

A h =p(S s nB {il} ) 

There are at least one and at most n = #1 non empty atoms of generation 1. Besides 
p{Ss) = {] %ieI A ll . 

Given (ii,^) £ ^ x 7 we call atom of generation 2 to the following set: Ai x j 2 = p(A^ n 
B {i2})- 

There are at least one and at most n 2 = {if I) 2 non empty atoms of generation 2. Besides 
p 2 (Ss) = U^^e/ 2 

By induction, if defined the atoms {-4i 1 ,t 2 ,...,ip}(ii,...,i p )eZ*' 
of generation p, we define the atoms of generation p + 1; 

^-ii,i2,— ,ip,ip+i = P\-^ii,h,—,ip ^ 
Then: pP(S s ) = U(i 1 ,i a ,.„ I i p )6/» ^l.fe.-A.- 

There are at least one and at most n p = {ffl) p non empty atoms of generation p. Fix p > 1. 
Denote A p to the finite collection of atoms of generation p: A p = {-Ai 1 ,i 2 ,...,ip}(i 1 ,...,ip)e/3'- 

Remark: It is easy to check from the Definition above that for all p > 1, if A & A p +\ 
then there exists A' £ A p such that A C A' . Precisely A = Ai lj i 2 ___ ) i +1 C Ai 2) ,„ > ^ >+1 = A' 
We define the diameter d p of the atomization of generation p: d p = max^g^ diam(A) 
where diam(A) = sup^^/ \\V — W\\ denotes the diameter of the atom A and it is if 
A = 0. We observe that d p is not the diameter of the union of the atoms, even if they 
intersect. 

Lemma 11 Let < A < 1 the uniform contraction rate of Theorem HI There exists a 
constant K such that for any p > 1 the diameter d p of the atomization of generation p 
verifies: 

0<d p < KXP- 1 , lim p ^ +00 d p = 0. 
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Proof: Take K = dianxB. Let us prove the thesis by induction on p > 1. 
An atom A G A\ of generation 1 is the image by p of a set in B. Then diam(A) < K. 
Suppose that all the atoms of generation p have diameter smaller than KX P . Let us 
prove the thesis for p + 1. An atom A G A p , p > 1 is the image by p of some set in B, so 
^4 C p(B). Take an atom A' G «4p+i- By definition: A' = p(A n Bm) for some ^4 € «4 P 
and some i G I. 

As j4 n £>{j} C p(-B) R due to Theorem H] p is contractive there: 

diam(A') = supyi y2 &A , \\~V~ 1 -V 2 \\ = 
= snpui tU 2 eAnB{ . } \\p{U l ) - p{U 2 )\\ < 
< supui^^B^ X WU 1 - U 2 \\ = 
= Adiam(A n B {i} ) < Adiam(A) < 
< KXiXP- 1 ) = K X p . □ 

Recall the distance dist(V, A) of a point V E B to a non empty set A C B: dist(V, A) = 
sup w&A \\V — W\\. Denote A and dA to the closure and frontier respectively, of the set 
A. A basic classic result from Topology asserts that, in any compact connected metric 
space, the following properties hold, if V $ A: dist(V,A) = dist(V,A) = dist(V, dA) = 
maxw^QA \\V — W\\ = \\V — Uq\\ for some Uq G dA. 

Lemma 12 IfV G S s then dist(V,W) > 5/2. 

Proof: Due to the metric properties in compact connected spaces recalled above, it is 
enough to show that dist(V, 8V) > 5/2. By contradiction suppose that miu^Q-p \\V — U\\ < 
5/2. Then there exists U G dV such that: ||V - U\\ < 5/2. 

V G Bg and by definition of uniform stability, if \\V — W\\ < 5 then: 

\\p k (V)- p k {W)\\ <a Vk > 1 (16) 

In particular V) — p(f)|| < a. 

On the other hand U G dV. Applying Theorem [7] there exists W G B such that 

\\U- W\\ < 6/2, \\p(U) - p(W)\\ > 3a. 
Joining the inequalities above with the triangular property: ||V — W\\ < ||V — J7|| + ||Z7 — 
W\\<5, 

3a<\\p(U)-p(W)\\ < 
<\\p(U)-p(V)\\ + \\p(V)-p(W)\\< 
<a + \\p(V)-p(W)\\, 

\\V-W\\<5, 2a < \\p(V) -p{W)\\. These last two inequalities contradict jTB}. □ 

Lemma 13 IfV€ S$i D B^y, then for all W G B such that \\W — V|| < 5/2 it is verified 
W G />', , [. 

Proof: Consider the compact set K = B \ B^y. Then: 

U G K O U£B {{ y. 
We have V G" K. Call d = dist(V,.ff) = inf(/^B {l} — U\\. As recalled when defining the 
distance of a point to a set, there exists Uq G dK = dUuy C dV such that 
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d=\\V -U \\= distQ^dK) > dist(V, dV). 
Due to LemmafT^ldistlV. dV) > 5/2. In resume we have proved that: d > 5/2, \rdjj^B {i} \\V- 
U\\ > 5/2 Therefore, if ||V - W\\ < 5/2 then W € B {i} . □ 

Definition. Indivisibility of the atoms. An atom A £ A p is indivisible if its closure 
A C for some i £ I. From the definition of atom, p(A) is a unique atom of generation 
p + 1, i.e. ^4 does not divide when applying p. Besides p is continuous in B^y so, it is 
continuous in A. Besides, A € p(B), so applying Theorem 

If an atom A is indivisible then p is continuous and uniformly contractive in A. 
Proof of Theorem [9] 

By Lemma [11] there exists a generation p > 1 of atoms such that the diameter d p < 5/2. 
Applying Lemma [13] and recalling that the diameter of a set is the same diameter of its 
closure, each atom A 6 A p is indivisible. From the remark at the end of Definition 1101 the 
atoms of generation p + 1 and later are contained in the atoms of generation p. Then all 
of them are also indivisible. 

Fix some non empty atom A = A 1 € A p . As it is indivisible: ^ piA 1 ) € «4 p +i- But 
any atom of generation p + 1 is contained in an atom of generation p, so there exists 
A 2 6 Ap such that p{A l ) C A 2 . From the indivisibility, p is continuous in A , thus: 
p(A ) C A . Applying the same argument to A2 instead of A x there exists A 3 6 A p 
such that p(A 2 ) C A 3 . We deduce p 2 {A l ) C p(A 2 ) C A 3 where A 1 , A 2 , A 3 G ^4 P . The 
family A p is finite, so there is some first pair of integer numbers ko > 1,to > 1 such that 
^fc _ ^ro+fc f , g ^ We conclude that some subfamily of _4 P is related in a finite chain: 

tf)c3 2 ,.., P (l')C^ 1 r" ) 

...p(A kQ )cA k0+1 ,..., 
...piA^-^cA^ =A ko (17) 

All these atoms are indivisible by construction, so for each A 3 , j = 1, 2 . . . , < ro — 1 in 
the finite chain (I17D . the Poincare map p : A 5 1 1— )• is continuous and contractive, with 
a uniform contraction rate < A < 1. In resume we have p r ° : A 1— > A and 
\\p ro (V) - p ro (W)\\ < X ro \\V - W\\ for all V,W £ ~A k ° . 

The Banach Fixed Point Theorem (Lages Lima [1970]), states that in any compact metric 
space M, given a uniformly contractive map / such that f{M) C M, there exists and is 
unique a point x\ € M fixed by /: f(x±) = x±. Besides all orbits by future iterates of / 
have limit {xi}, i.e. limfc_j. +00 f k {x) = x\, Vx E M. We conclude that p ro in the compact 
metric space A k ° has a fixed point V 1 such that: 

V 1 £A^, p r °(V 1 ) = V 1 , 

lim k ^ +QO p kro (V) = Vi \/V € A k ° 
Therefore V 1 is periodic by p of period ro, and its orbit L\ = {V 1 , pfV 1 ), . . . , /3 r ° _1 y 1 } is 
the omega limit of all the points in {jj~ ko+r ° A 3 , in particular of those in A . 
From the definition of atom and from p{S$) C Sg each atom is contained in S$. By Lemma 
[13] and using that p is contractive in each of its continuity pieces intersected with p(B), we 
deduce that all the points W such that dist(p(W),\J k °^° +l A j ) < 5/2 have u{W) = L x . 
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Such points W form an open set M D l/J=to° _1 A Besides L i c UjLto° _1 A ^ Then ^ 
is a neighborhood of Li. Therefore the basin of attraction B(L\) of the periodic orbit L\ 
contains a neighborhood of itself, verifying the definition of limit cycle. 
In particular, B(L\) D A £ A p . The construction above can be done choosing any first 
atom A\ e A p . We deduce that for each A € A p there exists a limit cycle L C B' such 
that IcB(L). 

Define C = {L C B : L is a limit cycle, ^4 C -B(L) for some A £ A p }. It is a finite and 
not empty collection because the family A p of the atoms of generation p is finite and not 
empty. By construction the union of the basin of attractions of all the limit cycles in £ 
contains all the atoms of generation p. Therefore: 

pP(B s ) = \J A eA P A c U L eC B ( L )- For a11 U the limit set u i u ) = bJ (p Pi J J ))- Then B S C 
[J LeC B(L). □ " 

5. CONCLUSIONS 

We described and analyzed a general mathematical model of a network of n inhibitory 
pacemaker neurons interacting by synapsis without delay. The size of the network must 
be finite but as large as wanted. We found discontinuities in the dynamical system due 
to the synaptic coupling and proved that they imply the existence of chaotic orbits. The 
Poincare return map p to the (n— l)-dimensional section B in the phase space, results from 
considering the state of the system immediately after each spike. We proved topological 
and measure properties of p as mathematical tools to obtain the dynamical results. We 
classified the systems according to the measure of the chaotic set and proved that, even 
being this set never empty, if it has not full measure there exist stable points and thus, 
limit cycles. It is unknown if there exist C 1 systems of this model exhibiting a set of chaotic 
orbits with positive measure, although under some other additional hypothesis, there are 
known results about the genericity of the systems exhibiting only limit cycles. 
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